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ABSTRACT 

The implications of extra sterile neutrinos for the Double Chooz experiment is expounded. The 
so-called "3+2" mass spectrum with 2 sterile neutrinos mixed with the active ones, is still allowed 
by the global neutrino data including MiniBooNE. We probe its impact on the resultant reactor 
antineutrino signal at the near and far detector of the Double Chooz experiment. The oscilla- 
tions driven by the additional mass squared difference due to the sterile states bring an energy 
independent constant suppression at both the near and far detectors. We study to what extent 
the measurement of #13 would get affected due to the presence of sterile mixing. We also give the 
projected sensitivity that Double Chooz will have to constrain the extra mixing angles associated 
with the sterile states. 
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1 Introduction 



Determining the mixing angle #13 is the next priority in the field of neutrino oscillation physics. 
Discovery of a non-zero value for this mixing angle is a prerequisite for determining the two other 
ingredients of the neutrino mass matrix 1 , viz, the CP phase 5c p and the sign of the atmospheric 
neutrino mass squared difference sgn(/S.m\ 1 ) . Various experimental proposals have been put for- 
ward to resolve this perplexing issue and measure this hitherto unknown mixing angle. Accelerator 
based experiments such as T2K and NOuA [1] involving conventional neutrino beams from pion 
decays are expected to come up in the near future. Superbeam upgrades of these facilities are 
also being envisaged. Pure v e and/or v e fluxes from beta decay of highly accelerated radioactive 
ions stored in rings is called Betabeam and under ingenious experimental set-ups can give ex- 
tremely good sensitivity to 9 13 [2]. The ultimate neutrino oscillation machine of course would be 
the neutrino factory [3], which could provide unprecedented #13 sensitivity, if indeed this mixing 
angle turns out to be extremely small. The main neutrino oscillation channel probed in these 
experiments are either the — > v e conversion channel P^ e or the v e — > conversion channel 
P eM , which is popularly known as the "Golden Channel" . The conversion channels suffer from the 
intrinsic problem of parameter degeneracies, whereby one ends up with multiple fake solutions in 
addition to the true one. The three kinds of parameters degeneracies are the so-called (#13, 5cp) 
intrinsic degeneracy [4], the (sgn(Am 31 ),5cp) degeneracy [5], and the (#23, tt/2 — #23) degeneracy 
[6], and they together lead to a total eight-fold degeneracy [7] of parameter values. 

The oscillation channel completely free of parameter degeneracies is the v e — * u e survival 
channel, P ee . This channel can be effectively probed in reactor based experiments. The crucial issue 
which needs to be addressed for maximum # 13 sensitivity is reducing the systematic uncertainties. 
It has been widely accepted that the best way of achieving this is by performing the experiment 
with two (or more) detector set-up. In these experiments a "near" detector is placed very close to 
the reactor cores and another one farther away, at the baseline optimal for observing near-maximal 
oscillations driven by Am| x [10]. Double Chooz [11, 12] is one such up-coming experiment which 
proposes to use the Chooz-B nuclear power plant in France. The other proposals that are being 
considered include ANGRA in Brazil [13], Daya Bay in China [14], and RENO in South Korea 
[15]. The first idea for using a near-far detector set-up in reactor experiments for measuring # 13 
was put forth by the KR2DET collaboration in Russia [16]. However, this experiment was shelved 
due to various reasons. Two other very good proposals which were eventually turned down by 
funding agencies include the experimental proposal to use the Braidwood reactor in the U.S.A. 
[17] and KASKA in Japan [18]. 

There are two-fold advantage of the experimental set-up with reactor neutrinos. Firstly as 
mentioned above, these experiments use the P ee channel, which is free from the problems of 
parameter degeneracies. Secondly, there are no problems of matter effects [19, 20, 21] in these 
experiments. However, there are at least two other experimental scenarios where very large matter 
effects in P ee can be exploited to provide stringent tests of $13 as well as sgr^Am^). It has been 
realized that extremely large matter effects inside the supernova leaves an imprint on the resultant 

-"^The remaining unknowns, which include the absolute neutrino mass scale and the Majorana phases (if neutrinos 
are indeed Majorana particles) cannot be ascertained in neutrino oscillation experiments. 
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neutrino spectrum and hence in principle its possible to determine #13 and sg^Am^) from the 
neutrino signal of a future galactic supernova [8, 9]. Very recently there has been a suggestion 
that in very long baseline experiments involving pure v e /v e fluxes such as Betabeams, one could 
effectively use the very large matter effects in P ee to pin down # 13 and sgn(Am\]) [22]. 

The recent declaration of the MiniBooNE results [23] has brought back the issue of sterile 
neutrinos on the front line of neutrino physics. The MiniBooNE experiment was designed to 
test the oscillation claim of the LSND experiment [24]. Since the oscillation interpretation of the 
LSND signal demands a mass squared difference Am 2 ~ eV 2 , it cannot be accommodated along 
with the solar and atmospheric data within a three flavor framework. One therefore needs one or 
more extra neutrinos and these species necessarily must be sterile. For one extra sterile neutrino 
one would get the so-called 2+2 and 3+1 schemes [25]. Even before the MiniBooNE results, the 
2+2 scheme was already disfavored from the solar and atmospheric data, while the 3+1 scheme 
suffered from severe tension between the LSND signal and the null results observed at other short 
baseline experiments [26]. However, all data could be fitted if one allows for 2 sterile neutrinos 
mixed with the active ones, leading to the so-called 3+2 mass schemes [9, 27, 28]. MiniBooNE has 
reported null signal in the energy range where it would have expected to see an excess of electron 
events if LSND experiment were due to two flavor oscillations. Though on the face of it this may 
look like a strong signal against the hypothesis of sterile neutrinos, in reality there is still ample 
room for their existence. For neutrino mass spectra with sterile neutrinos which predict Am 2 ~ 
eV 2 , one can see from the detailed analysis presented in [29] that even though the tension between 
the positive signal at LSND and negative signal at all other short baseline experiment including 
MiniBooNE, rules out the 3+1 scheme, the 3+2 scheme with 2 sterile neutrinos is still viable if 
one allows for CP violating phases. One should also keep in mind that even if it was proved that 
LSND was indeed wrong, it would only mean that sterile neutrinos in the eV 2 mass regime are 
ruled out. One could still envisage heavier sterile neutrinos which would have implications for 
astrophysics and cosmology, and these neutrinos could be mixed with active neutrino species. 

In this paper we look at the implications of these extra sterile neutrinos for the Double Chooz 
experiment. For concreteness we work within the allowed framework of the 3+2 scheme which gives 
a viable explanation of current world neutrino data. However, our results can be easily extended 
to other scenarios since the oscillation driven by the extra Am 2 corresponding to the sterile states 
anyway average out for the Double Chooz experiment. Our results can also be easily extended to 
the other reactor experiments mentioned before. We begin by discussing the oscillation probability 
in the 3+2 scheme in section 2. In section 3 we present the expected events in the near and far 
detector of Double Chooz if the 3+2 scheme was correct. Section 4 has our results with the full 
statistical analysis of the projected data set of Double Chooz. We end in section 5 with discussions 
and conclusions. 
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Figure 1: The v e survival probability as a function of energy for the near (left panel) and far (right 
panel) detectors of Double Chooz. The 6 different line types are for 6 combinations of sin 2 # 13 and 
sin 2 #i4 which are shown in the panels. We assume that sin 2 #15 = sin 2 #14. 



2 Oscillation Probability in the 3+2 Scheme 

The probability channel relevant for Double Chooz is P g g which for the 3+2 mass spectrum is 
given by 2 

P« = 1 -4El^| 2 |^| 2 -n 2 (^) -4El^| 2 |^| 2 sin 2 (^) , (1) 

i>j \ I k>l \ I 

where II s is the mixing matrix and the indices i and j run from 1-3, I runs from 1-5 and k s could 
be either 4 or 5. The oscillatory part of the third term of course will average out to 1/2 for Double 
Chooz, reducing the probability to 

P- e - e = 1 - 4£ mm/ sin 2 ~ 2£ TOW ■ (2) 

For the 3+2 mass spectrum we would have a 5 x 5 mixing matrix for which we choose the convention 

u s = R(e 45 )R(6 3 ,)R(e 34 )R(e 25 )R(e 24 )R(e 15 )R(e u )R(e 23 )R(e 13 )R(e 12 ) , (3) 

2 Note that for the Double Chooz experiment there are neither any matter effects nor any CP violating effect 
due the phases. The survival probability is therefore the same for neutrinos and antineutrinos. 
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where -R(%) are the rotation matrices and 9ij the mixing angle. We do not show the CP phases 
in Eq. (3) for simplicity. The mixing matrix with the above convention is expressed as 

/C15C14C13C12 C15C14C13S12 C15C14S13 C15S14 Si 5 \ 



U s = 



V 



(4) 



/ 



where Cij = cos%, Sij = sin%, and we show explicitly only the first row since the probability P ss , 
involves only them. We can see that apart from the two usual mixing angles 9 12 and #13 which 
appear in P ss with standard three generation oscillations, we have 2 additional angles, #14 and 
9 15 , which will affect the probability. If these mixing angles were zero, one would get back the P g g 
predicted by three generation oscillations. Of course for the energy and baseline of Double Chooz, 
the oscillations due to Am^ are extremely weak and as a result so is the dependence on 9 U . 
Therefore the dominant dependence of the probability would be on the mass squared difference 
Amj! and the mixing angles #13, 9 14 and #15. 

Current 3a constraints [30] on the parameters driving the leading oscillations in solar [31], 
atmospheric [32], K2K [33], MINOS [34] and KamLAND [35] are 



7.2 x l(T 5 eV 2 < Am 2 21 < 9.2 x 10" 



(5) 



0.25 < sin 2 9 X2 < 0.39 , 



(6) 



2.0 x 10~ 3 eV 2 < At/4 < 3.2 x 10~ 3 eV 2 



(7) 



sin 2 29 23 > 0.9 . 



(8) 



The best limit on the mixing angle 9 13 comes from the combined constraints from global oscillation 
data including CHOOZ [36] and is given as [30] 



sin 2 #13 < 0.044 . 



(9) 



The sterile sector receives constraints from the short baseline reactor and accelerator based experi- 
ments [26] which reported null signal, the LSND experiment [24] and MiniBooNE [23]. The best-fit 
values for the mass squared difference Ara^ and Am^ are 0.87 eV 2 and 1.91 eV 2 respectively, if 
the low energy MiniBooNE data is also included [29]. The best-fit values for the elements U e 4 and 
U e5 of the mixing matrix are 0.12 and 0.11 respectively [29]. Since in the convention adopted in 
this paper U e5 = sin 6*15 and f/ e4 = cos #15 sin 14 , this would translate to the best-fit values for the 
mixing angles as sin 2 9±4 = 0.012 and sin 2 #15 = 0.015. 

The Double Chooz experiment is being built to probe the mixing angle #13. However, since 
the probability depends also on the sterile mixing angles 9u and #15, these angles also can be 
constrained in this experiment. The parameters Am^ and Am^ are of course averaged out and 
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hence cannot be probed, and as discussed before, the solar parameters Am^ and sin 2 6*12 bring 
in a weak effect on P EE . In what follows, we will keep Am| x and sin 2 # 12 fixed at their best-fit 
values given in Eqs. (5) and (6). The true value of Am 2 x will be assumed to be 2.5 x 1CT 3 eV 2 
throughout. Also, just for the sake of simplicity we will take sin 2 9 U = sin 2 9 15 everywhere. 

In Fig. I we show the probability P g g as a function of the antineutrino energy for 6 different 
choices of the mixing angles #13 and 9 14. For all other oscillation parameters we stick to the 
assumptions mentioned in the previous paragraph. The left panel of the figure shows P g g at a 
distance of 280 m, which will be the average distance of the near detector from the 2 cores of the 
Chooz-B reactor. The right panel shows P gg for the far detector which will be at an average distance 
of 1.05 km from the 2 reactor cores. For the near detector, when # 14 = = #15, we expect P ss to 
be almost 1 at high energies. However at lower energies, for E ^ 4 MeV there is a Am| x driven 
dip in the survival probability, the extent of the dip depending on the value of sin 2 26* 13 . When 
we put #i3 = and allow 6*14 and 9 15 to be non-zero, we get an energy independent suppression, 
depending on the value of these mixing angles. When all three mixing angles are non-zero we have 
an energy independent suppression driven by 9 14 and 9 15 superposed on the energy dependent dip 
at low energy due to #13. At the far detector the Am^ driven oscillations of course are absolutely 
pronounced and the amplitude of the oscillations are determined by the value of sin 2 29 13 . The 
energy independent average oscillation due to the extra large mass squared differences due to the 
sterile states are superimposed on the standard oscillations. We can see from the figure that more 
than one combination of #13 and 9 14 (and #15) would give the same total suppression of the flux 
due to oscillations at the far detector. However, the shape of the resultant flux at the detector is 
expected to be different for the #13 and 9 14 (and 15 ) dependent oscillations. 



3 Results 



3.1 Number of Events at the Near and Far Detector 

The Chooz-B reactor complex consists of two reactor cores, east (Re) and west (%), with a total 
thermal power of 4.27 GW each. The v e flux at the near (N) and far (F) detector can be written 

as 

*(£) = T^£Af (10) 

where L N , F is the distance from the reactor core to the near (N) or far (F) detector 3 , Np s are the 
number of fissions per second for the isotope % in the reactor which we take from [37] and <fii(E) 
gives the corresponding energy spectrum 

<j) z (E) = exp(j2a kl E k ) . (11) 

\fe=0 / 

3 The far detector is at a distance 1114.6 m and 997.9 m from Re and Rw respectively, while the near detector 
is 290.7 m and 260.3 m away from them. 
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Figure 2: Number of events in 5 years expected at the near (left panel) and far (right panel) 
detectors of Double Chooz, as a function of the "visible energy" E vis of the detected positron. 
The 5 different line types are for 5 combinations of sin 2 #13 and sin 2 # 14 which are shown in the 
panels. We assume that sin 2 6*15 = sin 2 #14. 
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Figure 3: Total number of events in 5 years expected at the near (left panel) and far (right panel) 
detectors of Double Chooz, as a function of sin 2 # 13 . The black solid line is for sin 2 # 14 = 0, red 
dashed for sin 2 14 = 0.01 and green dot-dashed for sin 2 6^ = 0.02. We assume that sin 2 6^5 = 
sin 2 6*i 4 . 
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We assume a second order polynomial parameterization of a^i for the four isotopes (i = 235 U, 
239 Pu, 241 Pu, 238 U). The coefficients are taken from Table 2 of [38]. The v e are detected through 
their capture on protons 

u e +p^e + + n (12) 

whereby the e + and n give the prompt and delayed signal respectively in coincidence. The near 
and far detectors will be almost identical, consisting of a target volume of 10.32 m 3 of liquid 
scintillator comprising of 80% dodecane and 20% PXE, leading to 6.79 x 10 29 free target protons. 
The scintillator is doped with 0.1% gadolinium. The computed number of positron events in n th 
energy bin in the detector is given by 

N^ F = F R / dE ms / dE cr{E) R(E, E ms ) P SS (E, L NF ) £ ! ffi ' , (13) 

JE n Jo i 47TL^ F 

where E V i S is the measured visible energy of the emitted positron when the true visible energy 
El is ~ E - 0.8 MeV, with E being the energy of the incoming reactor antineutrino. The reaction 
cross-section is denoted by cr(E) and R(E, E vis ) is the energy resolution function of the detector. 
The quantity F R is given by 

F R = G x P x N p x T x (1 - d D ) x e D , (14) 

where, G is the Global Load factor (reactor efficiency), P is the Reactor thermal power, N p are 
the number of protons in the target volume, T is the exposure time, do is the dead time fraction 
of Detector D and en is the detector efficiency. Values of all quantities needed for calculating the 
positron rate in the detector are taken from [12]. 

We show in Fig. 2 the number of events expected in 5 years as a function of the visible positron 
energy, for the near (left panel) and far detector (right panel). We show the event spectra for 5 
different combinations of 9i 3 and 9 U (and #15) values. For the near detector it is mainly 9 14 (and 
#15) values which cause difference to the event spectra. At the far detector oscillations with both 
frequencies are important and hence all the three mixing angles make a difference. We stress that 
even though we do not show the statistical errorbars on this figure for the sake of clarity, one 
can easily check that most of the cases of mixing parameters displayed on this figure should be 
statistically distinguishable by combining the near and far event spectra at Double Chooz. 

In Fig. 3 we show the total observed positron events in 5 years as a function of sin 2 #13, for 
the near (left panel) and far detector (right panel). The three different line types in either of the 
panels show the results for a different choice for the value of 9 U (and #15). The net suppression 
of course increases as sin 2 9 U increases. This is true at both the near and far detectors. However, 
the sin 2 6*13 dependence of the suppression is extremely mild at the near detector and very large 
at the far detector. As expected, the largest suppression comes at sin 2 6>i 3 = 0.5 at which we have 
maximal oscillations. The most important thing we can note from this figure is that for a certain 
range of events at the far detector, any given observed rate could be predicted by a wide set of 
possible values of 6*13 and 9u (and #15). All these would then constitute degenerate solutions. For 
instance, we can see from the figure that if the far detector was to observed 7 x 10 4 events in 5 
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Figure 4: The 90%, 95%, 99% and 99.73% C.L. contours in the sin 2 13 - sin 2 14 plane for different 
assumed true values of the mixing angles marked in the figure by stars. Left panel is for the data 
from near detector only while the right panel is for data from far detector alone. We assume that 
sin 2 6*i5 = sin 2 6*14. 



years, then this would allow the sets of (0 13 , M ) values, (0.043, 0.00), (0.026, 0.01) and (0.008, 
0.02), as possible solutions. These would be degenerate solutions in this case. However, we can see 
by looking at the left panel that each of these combinations would predict different total rates at 
the near detector. Therefore by combining the near and far detected event rates, one can overcome 
this degeneracy problem. 



3.2 The #13 sensitivity 

For estimating the projected sensitivity of the Double Chooz experiment to the mixing angles #13, 
014 and 0i5, we perform statistical analysis of 5 years prospective data at both the near and far 
detectors. We define a \ 2 function on the lines of [12] 

[(1 + a corr + a° ncorr + a* )JV£(e /a ) - Nh(e true )} 2 



2 ir\ \ LV-"- 1 ^corr > ^uncorr ' sped 

X R (®true) = 2^ 1^ 



D=N,F i 



N i D + B i D + (AVfcm) 2 + (BhVbkd) 2 



n 2 (n D \ 2 (a* ) 2 

i "'corr | \ V^uncorr) , \ " V sped ) /ir\ 

"I 2 l ~ 2-^i 7 \2 7 \2 ' v ' 

&corr D=N,F V^uncorr) j {Aspect) 

where D runs for the number of detectors (near and far) and i runs over the number of bins, Qtme 
are the set of oscillation parameters at which the data is generated and fn is the corresponding 
set in theory. The systematic errors taken into account can be broadly characterized as normal- 
ization errors and spectral shape errors. The normalization errors include the 2% overall reactor 
antineutrino flux uncertainty which is relevant for both the detectors and therefore correlated. We 
denote this by a corr . The reminiscent uncorrelated normalization error between the detectors is 
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Figure 5: The 90%, 95%, 99% and 99.73% C.L. contours in the sin 2 #i 3 — sin 2 6*14 plane for assumed 
true values of the mixing angles marked in the figure by stars. Data from both near and far 
detectors are combined. We assume that sin 2 #15 = sin 2 6*14. 
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Figure 6: Sensitivity plot showing the A% 2 as a function of sin 2 #13 for the far detector only 
(left panel) and when near and far data sets are combined (right panel). Data is generated at 
sin 2 9 13 = and sin 2 #14 = for all curves in both panels. In the fit we fix sin 2 #14 = (black solid 
line), 0.02 (red dashed line) and 0.05 (green dot-dashed line). We assume that sin 2 9 15 = sin 2 6*14. 
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Figure 7: Sensitivity plot showing the A% 2 as a function of sin 2 #13 for the far detector only 
(left panel) and when near and far data sets are combined (right panel). Data is generated at 
sin 2 6*i3 = and sin 2 9 U = (black solid line) 0.02 (red dashed line) and 0.05 (green dot-dashed 
line). In the fit we fix sin 2 6*14 = for all curves. We assume that sin 2 6*15 = sin 2 #14. 
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Figure 8: Sensitivity plot showing the A% 2 as a function of sin 2 #13 for the far detector only 
(left panel) and when near and far data sets are combined (right panel). Data is generated at 
sin 2 #13 = and sin 2 #14 = (black solid line) 0.02 (red dashed line) and 0.05 (green dot-dashed 
line) . In the fit sin 2 #14 is allowed to vary freely and take any possible value for all curves. We 
assume that sin 2 6*15 = sin 2 #14. 
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taken as 0.6% and is denoted as cr uncorr . The uncertainty in the reactor u e spectral flux is denoted 
by a spect and we take 2% as its estimated value. This error is totally uncorrelated between the 
near and far detectors. There is also a bin-to-bin uncorrelated systematic error which is denoted 
by (J hin an d taken as 1%. We also include the background subtraction error <r bkd which is taken as 
1% of the background, B l D , and we assume that there is 1% background. We have checked that 
the background and its corresponding error makes almost no difference to our final results. For 
each set of oscillation parameter value taken in the fit, the function Xij(©tr«e) is minimized with 
respect to the parameters a corr , a^ ncorr and a % spect1 which are allowed to vary freely. 

In Eq. (15) we have explicitly put the subscript R on the x 2 function to denote the contribution 
coming from the Double Chooz set-up alone. In our numerical analysis, we have also included a 
"prior" on the allowed values of Am^ since we expect that the uncertainty on this parameter will 
see some reduction by the time the Double Chooz results are declared. Our full x 2 is therefore 
given by 

X 2 = x\ + Xlrior , (16) 

where, 

Xprior \ )i V / 

where we assume that Am^itrue) = 2.5 x 10~ 3 eV 2 and a^ m 2^ is 10% of Am^itrue). The 
total Xtot i s th en minimized with respect to some or all the oscillation parameters, to obtain the 
best-fit values, sensitivity limits and C.L. contours. In all our results presented in this paper, we 
marginalize over Am^. We keep Am^ and sin 2 6> 12 fixed at their best-fit value due to reasons 
discussed before. For the mixing angles #13, 9 U and 6*15, we will always mention whether they are 
free or fixed in the fit. 

In Fig. 4 we show the expected C.L. contours in the sin 2 #13 — sin 2 #14 plane. The left panel 
shows contours expected from analysis of data from the near detector only, while right panel shows 
the corresponding results when only the far detector data is analyzed. The points at which the data 
were generated are shown by the "star" marks in the figure. The different color shades show the 
90%, 95%, 99% and 99.73% contours. The Ax 2 for the C.L. contours correspond to 2 parameters. 
We see that the far detector can simultaneously constrain #13 and the sterile mixing angles #14 and 
9 15 . If the true value of sin 2 # 13 = 0.02 and sin 2 6> 14 = 0.05, we could measure the sterile mixing 
angles within the range, 0.025 < sin 2 9 U < 0.072, at 3a. If the true value was sin 2 9 13 = 0.02 
and sin 2 9 14 = 0.01 the lower limit for sin 2 9 U would be restricted by from below so that we 
would have 0.0 < sin 2 14 < 0.03. For sin 2 9 13 we see that at 95% C.L. 4 0.0026 < sin 2 9 13 < 0.037 
when data is at sin 2 #13 = 0.02 and sin 2 #14 = 0.05, and 0.003 < sin 2 #13 < 0.036 when data is 
at sin 2 9 13 = 0.02 and sin 2 14 = 0.01. We note from this that the limits on sin 2 # 13 at the far 
detector depends on the value of sin 2 9 U , albeit very slightly. Indeed the tilt of the C.L. contours 
towards the left shows a mild anticorrelation between the two mixing angles. We remind the 
reader that we keep sin 2 615 = sin 2 9u fixed throughout this paper. The contours for the near 

4 We give the 95% C.L. for sin 2 13 since above this the contours cross the y-axis. 
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detector explicitly show that there is hardly any sensitivity to sin 2 9± 3 , unless its value was very 
large, which anyway is already disfavored. It could however restrict at 3cr the sterile mixing to 
0.03 < sin 2 14 < 0.069 if the true value was sin 2 9 14 = 0.05 and 0.0 < sin 2 9 14 < 0.028 if the true 
value was sin 2 9 14 = 0.01. We see that the sterile mixing angle can be determined pretty well in 
either the near or far detector, with the precision in near detector being better due to its larger 
statistics. 

In Fig. 5 we show the allowed areas at 90%, 95%, 99% and 99.73% C.L. for 2 parameter fit 
when data from near and far detector are combined using the full expression given by Eq. (15). For 
the case where the data was simulated at sin 2 9± 3 = 0.02 and sin 2 14 = 0.05, the 3<r limits on the 
mixing angles are 0.007 < sin 2 9 13 < 0.037 and 0.031 < sin 2 9 14 < 0.068. For data at sin 2 9 13 = 0.02 
and sin 2 9 14 = 0.01, corresponding limits are 0.008 < sin 2 9 13 < 0.036 and 0.0 < sin 2 9 14 < 0.028 
We stress that we have given here the 3cr limits for sin 2 9\ 3 while in the previous paragraph they 
were at the 95% C.L. Note that the sensitivity to sin 2 9i 3 improves dramatically. The reason of 
course is that the near-far detector combination effectively assures that the correlated systematic 
errors cancel out, leaving just the small uncorrelated part. The uncertainty on sin 2 9± 4 and sin 2 #15 
is also somewhat reduced. Another significant feature we notice from the figure is that the near- 
far detector combination breaks the anticorrelation between sin 2 6>i 3 and sin 2 9\ 4 seen for the far 
detector in Fig. 4. The reason for this was discussed in the previous section. 

If the true value of #13 turns out to be very small or even zero then we would not see a positive 
signal at Double Chooz, at least in the three generation picture. We could nonetheless use the 
data to put an upper limit on the value of 9\ 3 at a given C.L. In what follows, we will give the 
"sin 2 #13 sensitivity reach" for Double Chooz. For this we generate the data at sin 2 9\ 3 = and fit 
it back allowing for non-zero value for sin 2 6*13. The resultant sin 2 9\ 3 sensitivity reach is shown in 
Fig. 6, 7 and 8. 

In Fig. 6 we have generated the data assuming the true values of sin 2 9± 3 = and sin 2 9\ 4 = 
(and sin 2 9 15 = 0). We find the A% 2 = (x 2 (sin 2 # i3 ) — xLin) f° r ever y value of sin 2 6* 13 , keeping 
sin 2 9 14 (and sin 2 #15) fixed at (solid black line) 0.02 (red dashed line) and 0.05 (green dot-dashed 
line) in the fit. We note that the A% 2 increases sharply with sin 2 9 U for the far detector alone. 
However, once we combine the data sets from near and far detectors, the effect of sin $14 becomes 
very small. This is related to the fact discussed before that the combined near and far data sets 
are almost independent of the correlation between sin 2 #13 and sin 2 9\ 4 and can measure both of 
them independently. In other words, the overall suppression or normalization factor due to sin 2 9\ 4 
gets canceled out from the fit for sin 2 9 13 , when data from the near detector is added to the data 
from the far detector. 

In Fig. 7 we generate the data assuming the true value of sin 2 9\ 3 = and taking sin 2 9\ 4 = 
(black solid line), 0.02 (red dashed line) and 0.05 (green dot-dashed line). Here we find the A% 2 
for every value of sin 2 #i 3 , keeping sin 2 9u (and sin 2 #15) fixed at in the fit for all cases. This 
situation might easily arise in practice if nature has indeed chosen a non-zero value for the sterile 
mixing angles and we being ignorant of that, try to fit the data by assuming that there were 
no sterile neutrinos mixed with the active ones. For the case where true value of sin 2 9± 4 was 
indeed 0, we have the standard three generation case and here we recover the projected sin 2 9\ 3 
sensitivity for Double Chooz. If only data from far detector was taken then we could put the 
limit sin 2 0i3 < 0.0095 (sin 2 2#i 3 < 0.038) at 90% C.L. and sin 2 ^ < 0.017 (sin 2 29 13 < 0.047) at 
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3(7. For combined data from near and far detectors we could restrict the angle to sin 2 9 13 < 0.006 
(sin 2 20i3 < 0.023) at 90% C.L. and sin 2 13 < 0.011 (sin 2 29 13 < 0.043) at 3a. However, if the true 
value of sin 2 # 14 was non-zero, we would have some difference in the sensitivity limit, if we were 
using data from the far detector alone. We would get a signal at the detector due to the sterile 
neutrinos and might confuse it with a signal due to sin 2 #13. It would look like we have observed 
a non-zero value of sin 2 6*13. We note from Fig. 7 that the best-fit sin 2 9 13 could be as large as 
sin 2 # 13 = 0.012 (0.027) if the true value of sin 2 14 was 0.02 (0.05). This is the fake solution that 
we discussed about in the previous section. Once we combine the data sets from both near and far 
detectors, the impact of the sterile mixing angle is negated to a large extent. However, for very 
large values of sin 2 #14 such as 0.05, we find that some residual confusion and the fake solution 
remains. Also, the upper limit on sin 2 9 13 at a given C.L. turns out to be different compared to 
the case where sin 2 9u was indeed zero. 

Fig. 8 is similar to Fig. 7 in most respect, except that here we allow sin 2 14 (and sin 2 9 15 ) 
to take all possible values in the fit. This would be the most democratic approach whereby all 
possible mixing angles are accounted for in the fit and marginalized over, while putting restrictions 
on #i 3 . The effect of the sterile mixing angle on the allowed values of sin 2 9 13 in this case is not 
dramatic. This is mainly due to the fact that sin 2 9u and sin 2 6*15 can now take non-zero values in 
the fit to give lowest possible x 2 , while in Fig. 7 the only mixing angle we could fiddle with was 
sin 2 9 13 since sin 2 # 14 and sin 2 15 were fixed at 0. When data at only the far detector is taken, the 
projected sensitivity on sin 2 #13 gets worse for true non-zero sin 2 #14, compared to when sin 2 9 14 
was truly zero in Nature. But once data from the near detector is added, the problem is solved 
and the standard projected sensitivity on sin 2 9±3 is restored. 

4 Conclusions 

The recent MiniBooNE results, even though disagree with the LSND data, leave ample room for 
existence of extra sterile neutrinos. In particular, the 3+2 neutrino mass spectra with 2 extra 
sterile neutrinos mixed with the active ones can satisfy the world neutrino data if one allows for 
CP violation. In this paper we probed the implications of the 3+2 mass spectrum for the Double 
Chooz reactor experiment. We showed how the event spectra at the near and far detectors of the 
Double Chooz experiment change when we allow for sterile mixing. The oscillations driven by the 
extra sterile neutrinos would produce a constant suppression at both the near and far detectors 
of Double Chooz. This is in contrast to the Am| x driven and #13 dependent oscillations which 
are almost absent in the near detector and only register their signal at the far detector. In the 
far detector, the extra oscillations due to sterile mixing would be superimposed over the normal 
flavor oscillation signal. In particular, we established that this sensitivity of the far detector to 
both 9 13 and the sterile mixing angles leads to a correlation between the 2 completely different 
kinds of mixing angles. 

We defined a x 2 function for the statistical analysis of the data and presented our results. We 
showed that the presence of sterile mixing angles alter, albeit slightly, the precision with which 
#13 could be determined in Double Chooz. The sterile mixing angle can be determined pretty 
well in either the near or far detector, with the precision in near detector being better due to 
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its larger statistics. The combined near and far detector data sets would be able to determine 
both mixing angles, though the precision in 13 goes down as the true value of the sterile mixing 
angle increases. We also studied how the "sensitivity to sin 2 # 13 " changes in presence of sterile 
neutrinos. We defined the sensitivity reach as the maximum value of sin 2 # 13 which would be able 
to fit the data at the chosen C.L., when its true value is exactly zero. If the true value of sterile 
mixing is indeed non-zero and we kept them fixed at zero in our fit, we would get significantly 
different results on the upper limit of sin 2 6> 13 , even after combining the results from both near 
and far detectors. This problem would stay if we allowed the sterile mixing to vary freely in the 
fit, but analyzed results from the far detector only. We finally showed that the problem could be 
fully negated only by combining results from both detectors and allowing the sterile mixings to 
take all possible values in the fit. 

In conclusion, presence of sterile neutrinos would leave its imprint on the signal at both the 
near and far detectors of the Double Chooz experiment. If only the far detector signal was 
considered, there is a possibility of confusing these sterile neutrino driven oscillations with active 
ones. However, by taking results from both detectors and allowing for the possibility of sterile 
mixing angles, one can probe both #13 and the sterile mixing angles correctly at Double Chooz. 
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